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Abstract — Restrictions on memory storage impose ul- 
timate limitations on the dimensionality of differential 
games problems for which optimal strategies can be 
computed via direct solution of the associated Hamilton- 
Jacobi-Isaacs equations. It is of interest therefore to 
explore whether, for certain specially structured dif- 
ferential games of interest, it is possible to decompose 
the original problem into a family of simpler, lower 
dimensional, differential games. In this paper we exhibit 
a class of single pursuer-multiple evader games for which 
a reduction in complexity of this nature is possible. The 
target set is expressed as a union of smaller, sub-target, 
sets. The individual differential games in the family are 
obtained from the original problem by taking as target 
set an element in the family of sub-target sets, in place 
of the original target set. We can exploit geometric 
features of the dynamic constraints and constraints of 
the problems arising in this way to reformulate them 
as lower dimensional, simpler to solve, problems. We 
give conditions under which the value function of the 
original can be characterized as the lower envelope of 
the value functions for the simpler problems and how 
optimal strategies can be constructed from those for 
the simpler problems. The methodology is illustrated by 
several examples. 

I. Introduction 

Interest in Pursuit-Evasion differential games (PE 
games) involving several players dates back to the 
1960's. Progress in this field is documented in the 
classical differential games literature, which includes 
the books by Isaacs [11], Pontryagin [15], Friedman 
[8], Krasovskii and Subbotin [14]. Constructing op- 
timal strategies for each player, finding the value 
of the game, deriving optimality conditions for the 
trajectories and establishing conditions for solvability 
of the game are typical objectives. 

The case of several agents, although it can be 
considered as a special case of the general framework 
for the problem, has been addressed separately by 
many authors. An early significant contribution was 
that of Pshenichnii [16], who studied the problem 
with many pursuers and equal speeds and derived 
necessary and sufficient conditions for solutions to 
this problem. Ivanov and Ledyaev [12] subsequently 
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studied the optimal pursuit in higher dimensional state 
spaces, with several pursuers and geometrical con- 
straints. Through the study of an auxiliary problem, 
related to the interaction between one pursuer and the 
evader and using a Lyapunov function, they obtained 
sufficient conditions of optimality. Chodun [4] and 
more recently Ibragimov [10] used the same approach 
to solve a collection of one to one problems (one 
pursuer one evader) for problems with simple dynamic 
constraints. 

Differential games, and multi-agent PE games in 
particular, have found applications in a variety of 
fields, for example in mathematical economics, where 
a game has been constructed to model the relation 
between agents [13], [6], in robotics where typically 
the emphasis is real time solutions and efficient com- 
putational methods. [9], [17]. 

Our proposed approach is to decompose the original 
game into a family of simpler, lower-dimensional 
games. This is achieve by expressing the original set 
as a union of smaller, target subsets. The individual 
differential games in the family result from the re- 
placing the original target set by each of the target 
subsets. Special geometric features of the dynamic 
constraints and constraints can be exploited to reduce 
the complexity of the new problems generated in this 
way. Making use of verification techniques originally 
proposed by Isaacs [5], [2], properties of viscosity 
solutions and techniques of nonsmooth analysis [1], 
[3], we give we a lower-envelope characterization of 
the value function, and show how to construct optimal 
strategies for the original problem from those for the 
simpler problems. A similar decomposition technique 
was used in [7]. 



II. The Hamilton-Jacobi-Isaacs approach to 

PURSUIT-EVASION GAMES 



The state y of a dynamic system, partitioned as n- 
vector components y = (yi, ■ ■ ■ , y m ), is governed by 



the equations 

y'i(t) = -gx{y{t)) ai {t) + hi(y(t))b(t) + h(y(t)) 

y'mit) = -9m+i{y(t))a m {t) 

+h m (y(t))b(t) + l m (y(t)) . 

(1) 

in which gi(.),hi(.) : K n ->■ K, i = 1, . . . , m + 1 and 
: E™ — > R n are given functions. 

For each i we interpret the state component yi to 
be the relative position of the evader with respect to 
the i'th pursuer. The dimension of the state vector y 
is therefore N — n x m. The n x m vector a = 
(oi, . . . ,a m ) comprises the n-vector pursuer controls 
and b is the n-vector evader control. 

The pursuer and evader controls a and b take values 
in the sets 

A = B n (0, Pa )x...x B n (0, Pa ), B = B n (0, Pb ) 

(2) 

for some given numbers p a , pb > 0. Define 

A := { meas. functions a : [0, +oo) — > A} (3) 
B := { meas. functions b : [0, +oo) — > _B} (4) 

It is assumed that 

m + 1 are Lipschitz 



(H): <?i(.),/ii(.U(-)> » = l,...n 
continuous, and 

9i{x)p a -h i (x)p b -\l i (x)\>0, Vx € 



iiA r 



Note that the maximum allowable magnitude of 
the velocity of the i'th pursuer depends on the states 
of all the pursuers and the evader, via the functions 
<?i(-)> and U{.). We may assume, without loss of 
generality that the g/s and hi's are positive functions, 
since the constraint sets A and B are symmetric about 
the origin. 



min \y t \<r}. 

ie{l,2...ro} 



Take the target set to be 

T={( yi ,y 2 ,-:y m )eR N 

jfci i,z...m f 

(5) 

in which r > is a specified number. Denote by 
yx(-) = Vx(-, a ,b) the solution of ([TJ, for given initial 
state x and controls a(.) e ^4 and 6(.) G i3. Define the 
hitting time for (given x, a and b) to be 

min{i : y x (t; a, b) e T} 
+oo if y x {t;a,b) £ TV< 

The player with control a (comprising the pursuers) 
seeks to minimize the hitting time, while the evader 
player, with control b, seeks to maximize it. It is 
convenient to transform the hitting time cost by means 
of the mapping 



t x (a,b) 



(6) 



The cost becomes: 

J(x,a,b) = ijj{t x {a,b)) = [ e~ s ds. (8) 

The transformation modifies the value function, but 
leaves unaltered the optimal strategies, owing to the 
fact that the transformation is monotone. 

Full specification of the differential game requires 
precise description of the nature of the control strate- 
gies involved. For this purpose, we follow the El- 
liot/Kalton approach, based on the concept of 'non- 
anticipative' strategies. The set of such strategies for 
the a-player is 

T := {a : B -> A : t > 0, b(s) = 6(a) for all s < t 
implies a[b](s) = a[b](s) for all s < t} . (9) 

A := {13 : A -> B : t > 0, a(s) = a(s) for all s < t 
implies f3[a](s) = /3[a](s) for all s < t} . (10) 

The upper and lower values of the game are 

u + (x) :— sup inf J(x, a, (3[a\) 

u~ (x) := inf sup J(a[6], b) . 

Under the stated hypotheses, and in view of the 
fact the Isaac's condition is satisfied, the upper and 
lower values coincide, for arbitrary initial state x. The 
common value defines the value function u(.) for the 
game. Thus u(x) = u + (x) = u~(x) for all x. 

It can be shown that, under the state assumptions, 
the value function is the unique uniformly continuous 
function u{.) : M. N — > K, which vanishes on on T and 
which is a viscosity solution on R \T of the following 
Hamilton- Jacobi-Isaacs equation 

u(x) + H(x,Du(x)) = xeR N \T 



x G 

u(x) = x € dT 

where The Hamiltonian is 

H(x,p = (pi, . . . ,p m )) = 



(11) 



Pa 



i=l 



9i(x)\Pi\ ~ Pb\^K{x)Pi 



ip(u) 



1 — e u if u < +oo 
1 if u = +oo 



(7) 



i=i 

~y^ji(x) - Pi - 1 . 

Furthermore, the value function is locally Lipschitz 
continuous. 

Solving this equation for the value function of the 
game, yields the optimal strategy of each player for 
initial state xq as a(t) — S(y Xo (t)) and b(t) = 
w {yx (t)) where 

S{z) £ argmaxmin{(5(a;)a — h(x)b — l{x)) ■ Dv(x)} 



W(z) € argminrnax{(<?(iz;)a— h(x)b— l(x))- Dv(x)}. 



Example 1. We examine a simple, preliminary ex- 
ample. Consider the pursuit-evasion game with two 
pursuers p\ , p 2 and one evader e, and the positions of 
each evolve in ID space. We take the dynamics to be 



Pi 

P'-2 

e' 



§ai 



«2 
b 

Pi(o) 2 = P i 

Pa(0) 



(12) 



P°2 



e(0) 



where a u a 2 ,b € £(0, 1) = [-1, 1], Pi,pi, e £ E. The 
speed constraint on pursuers is assumed to be greater 
than that of the evader; this ensures that the optimal 
hitting time is finite for an arbitrary initial finite state. 

Let us now consider a reduced formulation. We 
translate the origin to the position of the e-player, 
so the z'th state variable becomes the position of the 
evader relative to that of the i'th pursuer. There results 
the reduced dynamics 



V 2 = -02 + I 

y 1 (0)=p o 1 -e° 
y 2 (0)=p° 2 -e° 



(13) 



here we have ai,a 2 , b G B(0 7 1), yi,y 2 G (-co, +oo). 
The target set T becomes the union of neighbourhoods 
the two axes: capture occurs when mirij e { 1 2 } \Pi — 
e\ < r, for some specified r > 0. 



pl 



p2 



al b a2 

Fig. 1. A monodimensional PE game with two pursuers. 



The HJI equation is 



Du(x)} 



u(x) + maxmin {(|&i — |, a 2 — \ 

\ xe [0,+oo] 2 \T 

u(x) =0 x G dl~ 

(14) 

Observe that, with this formulation, the dimension 
of the domain in which we seek to solve ( |T4] > is two. 
The rapid growth in complexity of the problem with 
increase in the number of the pursuers places severe 
restrictions on the computability of this equation for 
state dimensions greater than three. We now propose 
a decomposition technique to address this challenge. 



III. A DECOMPOSITION TECHNIQUE FOR THE TO 
PURSUER PROBLEM 

In this section we present a decomposition technique 
to overcome the complexity of the multiple pursuer 
game introduced in Section 2, for high state dimen- 
sions. The key idea is to decompose the target set T 
as a union of smaller sets %, i = 1, • • • , m: 

T t = {(yi,y 2 , -,y m ) g : < r}, i = 1, ...,m. 

(15) 

and to relate the value function u(.) for the original 
game to the value functions Ui(.) for the games in 
which the 77s replace T. The «i(.)'s are viscosity 
solutions to the HJI equation above, with modified 
boundary condition: 

Ui(x) + H(x,Dui(x)) =0 xeR N \% 
Ui{x) = x G Ti 

(16) 

Under the hypotheses, the Ui(.)'s are Lipschitz 
continuous. 

Define the index set 

I{x) := {i G {1, . . . , m} : Ui{x) = minu(i)} . 

i 

Theorem 3.1: Assume condition (C) is satisfied: 

(C): for arbitrary x G R N \T, any convex combi- 
nation {Xi | i G I(x)} and any collection of vectors 
{£i G d L Ui(x) | i G I(x)} we have 

Then 

u(x) = min{ui(x), . . . ,u m (x)} 

i 

for all x G R N \T. 

Remark 3.2: Here the limiting superdifferential 

d L Ui(x) is the set 

d L Ui(x) := lim sup d F Ui{x) , 

in which d F Ui(x) is the super Frechet differential 

d F Ui(x) 

■ TTT1./V 1- U i (x')—U i (x)—p-(x / — x) ^ n -| 

{p G R" : lim sup — \^- x \ 1 < 0} . 

x' — >x 

Condition (C) is automatically is satisfied if H(x, .) is 
a convex function, but is in fact significantly weaker. 

Proof Outline. Define 

u{x) :— min{uj(x); i G 1 . . . to} for all x (18) 

Our aim is to show that u(.) coincides with u(.). Since 
u(.) vanishes on T, and in view of the facts that u(.) 



is Lipschitz continuous and that u(.) is the unique 
uniformly continuous viscosity solution to the HJI 
equation satisfying this boundary condition, it suffices 
to show that u(.) is such a solution. 

That u(.) is a super (viscosity) solution follows 
directly from the definition of super solution and the 
fact that u(.) is the lower envelope of a finite number 
of super solutions. It remains to show then that u(.) is a 
sub solution. Take any x G M. N \T and any £ G d F (x). 
According to a well-known characterization of sub 
solutions, it suffices to demonstrate that 



u + H(x, f) < . 



(19) 



But, by the 'max rule' for limiting subdifferentials of 
Lipschitz continuous functions, applied to — u(x) — 
max(— u)(%), we know that, for some convex com- 

i 

bination {A^ | i £ I(x)} and set of vectors {£j G 
d L u l {x) \i e I(x)}. 

z = E x & ■ 

iei(x) 

For each i g I(x), there exists sequences Xj — > x and 
Cj — > £,i sucn that, f° r sach i, £ G d F (x l j) as j — » oo. 
Since is a subsolution, H(xp£,j) < 0. But then, 
since H(.,.) is continuous, 

<Um sup ff <0 

j 

It follows from condition ( fl7] > that 

H(x,0 = H(x,J2^i) < Y^*iH(x,& . 

i i 

then 

W + H X,J2 X iPi < E + E X i H ( X > Pi) 
\ i / i i 

We have confirmed {19) . □ 
We now discuss the role the decomposition tech- 



nique, summarized as Theorem 3.1 in reducing the 
complexity of the differential game. 

Consider again Example 1. The target can be ex- 
pressed as a union of two sets: T := 71 U T2, in which 
% := {(xx,x 2 ) € M 2 : \xi\ < r}. For the target 71 
(respectively T2), the value function is clearly inde- 
pendent of X2 (respectively x\). So g|-iii(xi, x%) = 
and gf-it 2 (xi, X2) = 0. m(.) therefore satisfies 



Ui + max mm 

ai b 



in < - 



i-hi 



+ 1) 



1 



Ml = 



3 1 * 1 1 21 9xi 

afi G (r,+oo],a;2 G [0, +00) 
sci G [0,r],x 2 G [0,+oo) 

(21) 



Fig. 2. Some level sets for the value function of the monodimen- 
sional problem with two pursuers 



where a\ G B n (0,p a ) and b G B n (0,pb). This is a 
ID equation for a fixed X2 and is constant for a fixed 
a?i. It has solution 

mix) = 1 - e - 6{xi - r \ ip^im) = 6(n - r) (22) 

Similarly 

u 2 (x) = 1 - e - 2 C*»-»0, V" 1 ^) = 2(x 2 - r) . (23) 

The gradients of these two functions at x are of the 
form (fci(x),0), (0,^2(2;)) for non-negative functions 
fci(.) and &2(.). It is simple to check that condition (C) 
is satisfied. According to Thm. |3.1| the value function 
as lower envelope of and U2(-)> thus 



— 6(xi — r) 
-2(x a -r) 



if xi < |x 2 + |r 
if xi > 7tx 2 + lr. 



(24) 

Remark 3.3: These optimal strategies are consistent 
with intuition. When one pursuer is very close to the 
evader, that pursuer's location alone affects the strategy 
of the evader, as expected. This feature of the solution 
is evident from the formulae for the value function, 
which reveal that, for states far from the X\ axis 
but close to the X2 (for example), the value function 
coincides with the function mix). 

The preceding analysis can be generalized to cover 
a unidimensional problem with m pursuers. 

Theorem 3.4: Assume hypotheses (H). Denote by 
: K — > K the solution of the following equation 

Vi{xi) +maxmin{/ i (a; l ,ai,6) • DviixA] = 1 

ai b 

X, G (r, +00] 
Xi G [0, r] 



v{xi) = 



where /^(xj, a, h b) := g l (x)a i - hi(xi)b - a* G 

B x (0,p o ), be BxiO, Pb ). 

Uiix = (xi,X2...Xi, ...X m )) = Wi(Xi) 



Then the value function u(.) is 

u{x) = mva{ux{x), . . . ,u m (x)} 

Proof: It is straightforward to confirm that Ui{x) 
is a viscosity solution of the problem ( [To} , using 
the fact that ffK^) = for a U * 7^ j- Then, in 
view of the uniqueness of viscosity solutions for ( fT6| ) 
([1] Theorem 3.1), we know that is the unique 
viscosity solution of the decomposed problem ( fT6| ). 

Making use of the fact that that G d L ui{x) has 
the structure (0, . . . , 0, fe,, 0, . . . , 0), where fc^ > 0, we 
deduce Then, for any G d L Ui(x) and £j G d L Uj{x) 
with i,j G /(x), we have 

+ £7) - 

Pa (ffiN +Pi|fcj|)-/96|/lifc i +/ljfc i |-i r (0...fe i ...fej...O) 
= (Pa5i| fc d - Pb|/l l fci| - k ■ (0...fej...0)) 
+ (Pa&N - Phl^'l - ' (0...fc 3 ...0)) 

= JJ(a:,&) + (25) 

this confirms condition (C) of Thm. |3.1| is satisfied. 
The stated characterization of the value function now 
follows from Thm. 13.11 

□ 

IV. Examples and numerical tests. 

In this section we solve some higher dimensional 
problems using the proposed decomposition technique. 
We note that memory storage constraints impose fun- 
damental limits on the state dimensions of problems 
for which solutions to the associated HJI equations can 
be computed directly. Indeed, MATLAB implemen- 
tations using a heap-based Java VM system are not 
feasible for N dimensional arrays, for N > 5 . This 
provides the motivation for studying decomposition 
techniques which are, in principle, applicable for very 
high dimensional problems. 

Test 1 

In this example we study pursuit-evasion game on a 
plane involving several pursuers and one evader, with 
the help of Thm. 3.1. 

Consider the problem in a reduced space where the 
variable Xi is the distance between the evader and the 
z'th pursuer and every pursuer has constant maximum 
speed gi(x) = 1. The evader has constant maximum 
speed hi(x) = 0.9. 

The i'th unidimensional problem is, for control 
constraint sets a. L G -Bi(0, 1) and b G Si (0,1), has 
associated HJI equation 

Vi{xi) + min max {((^a^ — hi(x)b) ■ DiVi(xi)} = 1 

a-i b 

< Xi G (r, +00] 

v(xi)=0 Xi&[0,r]. 

(26) 



We extend now the 1^ solution in the space of dimen- 
sion 5: Vi(x) : M 5 -> R is defined as 

u t (x = (xi, ...xi, ...x 5 )) = Vi(xi) . (27) 

The value function of the problem is simply 

u(x) = min{u!(a;), . . . ,u 5 (x)}. (28) 

Figs. [3] |4] we show simulations for two different 
with different starting points. 
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Fig. 3. Test 1: optimal trajectories at various times. A X indicates 
the point of capture. 

2r 
1.5 - 




-1.5 - 

-2" 1 1 1 1 1 1 1 1 1 >- 

-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



Fig. 4. Test 1: optimal trajectories at various times. A X indicates 
the point of capture. The optimal trajectory of the evader is waiting 
the capture. 

Test 2 




Fig. 5. Test 2: Target set for m = 2, n = 1, r = 0. 

In this example, we consider a PE game involving 
several pursuers and one evader, but now in presence of 
an obstacle (a river, a different kind of medium) which 
affects the velocities of the pursuers and the evader. In 
this case we cannot pass to the reduced coordinates, 
since this will result in a loss of the information about 
the position of the agents. Accordingly, the first m 
block components of the state are associated with the 
pursuers and the (m + l)'th block component with the 
evader. The positions of the agents, evolving in R n , 
are governed by the following equations 

' yiit) = -g(y(t)) ai (t) 
l/aW = -g(y(t))a 2 (t) 
< ••■ (29) 

y'mit) = -g(y(t))a m (t) 

. y' m+1 {t) = h{y{t))b{t), 

where the state y := (yf,y^,...,y^ l+1 ) T G and 
every j/ { e R n . We assume that g : R n — > R + and 
h : R n — > R + , so every pursuer has the same velocity 
rule (depending from the position; in the previous 
case y was the distance between agents). This will 
give us a geometric property which will guarantee the 
decomposability of the problem. The appropriate target 
set in the present setting is 

T={(yf,yt,...,yT n+1 )eR N : 

min \y t - y m+ i \ < r \ (30) 
»e{i,2...m} J 

Fig. [5] illustrates the target set in the case m = 2, 
n = 1. 

The Hamiltonian, which is 

m 

H{x,p) = g(x)^2 p a \pi\ - h(x)p b \p m+1 \ - 1 , 

i=l 



is not convex in the p variable. 

Theorem 4.1: Assume Q. Let u t (.) : R" -> K 
the solution of the following equation 

Vi{xi) + maxmin {(g(xi)cii, —h(x)b) ■ Dvi(xi)} = 1 

ai b 

x £ (r, +oo] 2 ™ 
v(xi) = Xi € {|afi - x 2 \ < r} 

(31) 

with ai € B n (0,p a ), b e B n (0,pb). Define the 
function Ui(x) : R N —> R as 

Ui(x= (xi,...,,x m )) =Vi(xi) . (32) 

Then the value function for the original problem is 

u(x) = mia{ui(x), . . .,u m (x)} 

Proof: To simplify the notation, suppose that n = 
1. (The case for a general n is treated in the same way). 
We can see that the Hamiltonian is convex in p along 
rays in every direction with exception of the e m +i 
direction. (Here, e, is i'th canonical basis vector). 

We can repeat the main steps in the proof of 
Theorem |3.4| with the exception of the verification of 
condition (C). In this case we know, for geometric 
reason, the superdifferential of the function u{.) is an 
element aligned with e m +i, i.e. 

Ci, € d F u(x), (& - • e m+1 = 0. 

This follows from the fact that e m +i is tangential to 
the switching interface. Writing Ui and Uj for two 
reduced value functions, we can show that itj(.) = 
it,-(.), i.e. they solve the same equation. It follows that 
the switching interface is located where two reduced 
value functions coincide, i.e. where Uj(x) = Ui(x). 
Writing n for the normal of the switching interface, 
we have 

n-e m+1 = ( 0, . . . , 0, 1" 1 , 0, . . . , 0, ~,0, ■ ■ ■ ,0 ) 
\ oxi axj J 

•(0,...,0,l) t = 0. (33) 

Condition (C) can now be validated, the state rep- 
resentation of u(.) is therefore valid by Thm. 3.1. 

□ 

As an example, consider the n — 2 case. An evader 
has position denoted by x e <E R 2 and m € N pursuers 
have positions denoted by x\, X2, ■■■X n £ R 2 . Take 
m = 3. It is assumed that 

J g(x) = 1 — 0.5cos(7rx) if < 0.5 
\ g(x) = 1 elsewhere 

and h(x) = 0.4. 

Figs. [6] |7] show some optimal trajectories and the 
level sets of the velocity function of the pursuers. 




Fig. 6. Test 2: optimal trjectories for a 3pursuers game. 



0.4 
0.2 



























( 






! E 










0< 




d t 

















1 1 ' 1 ' 1 1 1 ' ' L_ 

-1 -0.8 -0.6 -0.4 -0 2 0.2 0.4 0.6 0.8 



[4] W. Chodun, Differential games of evasion with many pursuers, 
/ Math. Anal. Appl. 142 (1989), no. 2, 370-389. 

[5] R.J. Elliott and NJ. Kalton. Values in differential games. Bull. 
Amer. Math. Soc. Vol. 78, No 3 (1972), 427-431. 

[6] G.M. Erickson, A differential game model of the marketing- 
operations interface, Eur. J. Oper. Res. 211 (2011) 294^402. 

[7] P. Falugi, C. Kountouriotis and R. Vinter, Differential Games 
Controllers That Confine a System to a Safe Region in the 
State Space, With Applications to Surge Tank Control, IEEE 
T. Automat. Contr. (2012) 57(11): 2778-2788 . 

[8] A. Friedman, Differential Games, John Wiley & Sons, New 
York, USA, 1971. 

[9] J.Hu and M. Wellman. Multiagent Reinforcement Learning: 
Theoretical Framework and an Algorithm. Proceedings of 
ICML, (1998) 242-250. 
[10] G.I. Ibragimov, Optimal pursuit of an evader by countably 

many pursuers. Differ. Equ. 41 (2005), no. 5, 627-635. 
[11] R. Isaacs, Differential Games, John Wiley & Sons, New York, 
USA, 1965. 

[12] R.R Ivanov and Yu. S. Ledyaev, Time optimality for the pursuit 

of several objects with simple motion in a differential game. 

Trudy Mat. Inst. Steklov. 158 (1981), 87-97. 
[13] S. J0rgensen, Optimal production, purchasing and pricing: A 

differential game approach, Eur. J. Oper. Res. 24 (1) (1986) 

64-76. 

[14] N.N. Krasovskii and A.I. Subbotin, Game-Theoretical Control 

Problems, Springer, New York, 1988. 
[15] L.S. Pontryagrin, Izbrange Trudy (Selected works), Moskow, 

Russia, 1988. 

[16] B.N. Pshenichnii, Simple pursuit by several objects, Cybern. 
Syst. Anal. 12 (3) (1976) 484^485. 

[17] R. Vidal, O. Shakernia, J. Kim, Associate Member, H. Shim, 
S.Sastry, Probabilistic Pursuit-Evasion Games: Theory, Im- 
plementation, and Experimental Evaluation, IEEE T. Robotic. 
Autom., (2002) Vol. 18, No. 5. 



Fig. 7. Test 2: optimal trjectories for a 3pursuers game. 
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